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^ ■ I. INTRODUCTION 

(N : 

' In the context of the "old" Standard Model (SM) with zero neutrino masses, the lepton flavor is conserved. As a 
. result, the LFV decays such as /i ^ eee and fi ~> are strictly forbidden within the "old" SM. Within the "new" SM 
rS i' with sources of LFV in the neutrino mass matrix, such decays are in principle allowed but their rates are suppressed by 
JL , powers of the neutrino mass and are therefore beyond the reach of any searches in the foreseeable future 1\. However, 
a variety of beyond SM scenarios can lead to rates for these processes exceeding the present experimental bounds Q : 

L—J, Br(/i ^ e7) < 1.2 X 10"" and Br(/^" ^ e"e+e") < 1.0 X 10"^2 . (1) 

Notice that the bound on fj, ^ eee is even stronger than the famous bound on /i 67. In the context of models 
. like R-parity conserving MSSM in which the new particles can appear only in even numbers in each vertex, these 
' processes can only take place at the loop level. Since /i — > eee is a three body decay, in such a model, Br(/i — > eee) is 
', suppressed relative to Br(/i ej) by a factor typically of order of e^/(167r^) log(m^/me) However, in the models 
fSj ' that new particles can appear in odd numbers at each vertex {e.g., in R-parity violating MSSM) the process /i —^ eee 
' can take place at tree level and as a result, its rate can even exceed that of /x — > e7 

' It is rather well-known that by measuring the angular distribution of the final particles relative to the spin of the 
initial muon in fi ej and a — > eee, one can derive information on the chiral nature of the effective Lagrangian 
. leading to this process In the case of /i ^ eee, as shown in the literature 0, HUli the angular distribution 

^ ' of the final particles relative to the spin of the initial muon also yields information on certain combinations of the 
• ^ , CP-violating phases. 

Recently, it has been shown in [^, [lol [Tlj that if we measure the polarization of the emitted particles in /i — > 67 and 
^ fiN — > eA^, we can derive information on the CP- violating parameters of the theory. It was pointed out in that 
by measuring the spin of the more energetic final positron in /i+ — > e"'"e~e"'", some information on the CP-violating 
phases can be derived. The analysis was performed in the framework of the models such as _R-parity conserving MSSM, 
in which the dominant contribution to fi eee comes from a penguin diagram (i.e., fj,~^ 7*e+ — > e~''e~e"*'). In this 
letter, we focus on the case that /x"*" e+e~e+ happens at the tree level through the exchange of heavy particles. We 
show that the transverse polarizations of the emitted particles in /i+ ^ e+e~e+ provide us with information on the 
CP-violating properties of the effective Lagrangian. The information derived from the spins of the emitted electrons 
and positrons are complementary to each other. In this letter, we focus on /z"*" — > e+e~e"'". Similar arguments hold 
for /i" e"e+e". 

We show that the method proposed here is sensitive to a combination of the phases in the effective Lagrangian 
which is different from those that can be derived by methods discussed in the literature [4j, 7, 8] . The effectiveness 
of each method depends on the relative magnitude of the different terms in the effective Lagrangian which in turn 
depends on the details of the underlying model. 
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II. THE RATE OF ^+ ^ 6+6 6+ 

Consider a general beyond SM scenario leading to /i+ e+e^e+. After integrating out the heavy states, the effect 
can be described by an effective Lagrangian of form 



where 



and 



i: = £i+£2 (2) 



r P ^- 1 + 75 ^^- 1 -75 X , p ^- 1-75 ^^- 1 + 75 ^ , 

^1 = e)(e — —e) + B2(ti — — e)(e ^ e) + 

^ ,_l + 75 w_l + 75 ^ , ^ ,_l-75 ^.-l-75 ^ 

+ Gi(Ai7''^^e)(e-7.i^e)+G2(A27''i-^e)(e-7.i^e) + H.c. (3) 



- ALfl[r,r]^^eF^. + A;^/i[7^7'^]i-^eF^. + H.C. (4) 



Notice that by using the identities {a^)af3{o'p,)^s = '^^a-y^ps and {a'^)ai3 = ^i3s{o'^)s-y^ia (where en = eoo = and 
£10 = — £01 = 1) and employing the fact that the fermions anti-commute, we can rewrite the terms on the first line of 
Eq. (O as 

Bl._l + 75 i.w_l-75 ^ ^2 1-75 .,_l + 75 ^ 

— 2"*-^ 2 ^ >s) - -ylA^ 2 ^ e)(e^— 7^e) . 

In the literature, it has been shown that by studying the angular distribution of the final particles relative to the 
spin of the initial muon, one can derive information on Re[Ai^i3* - AlB^], Re[^_RG'2 — AlGW, ItoIArBI + A^B^] 
and Im[ylflG'2 + AlG\] (see e.g., [1, 0, ^]). However, by this method the phases of Ci and C2 cannot be derived. 
Moreover, if and Ar are zero or suppressed, this method loses its effectiveness. 

By studying the energy spectrum of the emitted particles, it will be possible to differentiate between the different 
terms in £1 and £2- For example, the Al and A]^ couplings lead to a sharp peak in the distribution of the square of 
the invariant mass of a e~e''" pair [i.e., in the distribution of (-Pg- — Pe+Y] at (m^/2)^. Such a peak does not appear 
if the dominant contribution comes from Eq. The A^ and Ab, couplings are generally loop suppressed but, as 
we see below, terms in £1 can appear in a wide range of models at tree level. In the present letter, we only consider 
terms in £1. 

The Bi and Ci couplings in £1 can originate from the exchange of a heavy neutral scalar field (or fields). Let 
us demonstrate this through a simple toy model. Consider two heavy complex fields, 4>i and 02 with the following 
couplings 

-1 + 75 _1~75 -1 + 75 -1 — 75 

ff^iL^iM — ^ — e + 5^fl02Ai — ^ — e + geL^ie — - — e + ge_R<^2e — - — e + H.c. 

It is straightforward to show that the effective Lagrangian resulting from integrating out the heavy states (pi and 4)2 
is of form ([3]) with 

a\ — 5 — , 02 — 2 — ' <-^i — (-^2 — (J • 

ml m, 

(pl 02 

If we swap 01 and 02 in the third and fourth terms {i.e., taking £ = [5;iL0iA(l + 75)e + 5/jfi02M(l — 75)e + 
5eL02e(l + 75)e + gei?0ie(l - 75)e + H.c.]/2), we find Ci = g^L9tR/'ml^, C2 = g^flSeL/'^^s and Bi = B2 = 0. 
Taking 0i and 02 real, we find that Bi, B2, Ci and C2 are all nonzero. It can similarly be shown that Gi and G2 can 
originate from the exchange of a doubly charged scalar field. The R-parity violating MSSM at the tree level leads to 
couplings of form Bi and B2 That is while in this model, the and Ar couplings are loop suppressed. 

Notice that while the Bi couplings are chirality conserving, the Ci and Ci couplings are chirality-flipping. Notice 
that under the parity transformation, Bi and Ci transform as Bi <^ B2, Gi G2 and Gi 44> G2. Thus, \Bi ~ B2\, 
|Gi — G2I and |Gi — G2I can be considered as measures of the parity violation. Under the CP transformation 



Bi T]Bl , B2 => , Gi ^ ryGI , G2 ^ r/G^ , Gi ^ rjC*, and G2 ^ , 
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where ry is a pure phase that comes from the freedom in the definition of the CP-conjugate of the electron and the 
muon. By rephasing the electron and/or the muon field, either of these couplings can be made real. Thus, the 
Lagrangian in Eq. ^ contains five physical CP-violating phases. 

From the formulas derived in the appendix, we find that the differential decay rate of — > e^e^e^ is 



2v 



E 



dE, dE,+ dfl 



1287r4 



+ \B2f){-ml + 4mlE^+ + 3mlE, 



(|Bip-|B2nP^cos^?_(m4 
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2Elm^ - 
8ml E,E 
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(|C2p + 16|Gin(f 
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(5) 



where dQ is the differential solid angle determining the orientation of the emitted electron. We have summed over 
the spins of the final particles and integrated over 0+ which determines the azimuthal angle of the emitted positrons 
(see Fig. [J). 

Integrating over the energies of the final particles, we find 
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dT 



dcosil 



0.0208 . 
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1287r4 



: cos 0-{- 



i\Bi\'- 
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\Ci\ 
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8|Gi|2 



8(|G2p 



^8|G2p) 
IGiH) 



(6) 



The total rate of / {dT/dn)dn is given by 
Br(/z — > eee) (see Eq. [T]), we find 



\B2 



\Ci\ 



\C2\')/2) 



|Gip + IG2P). From the bound on 



\Bi\ 



IBo 



IC1P + IG2P 



K|Gi 



IG2P) 



< 



1 



(200 TeV)^ 



From Eq. we observe that the dependence of the coefficients of -I- |i32|^ and |Gi -|- IG2P -|- 16(|Gip -I- IG2 1^) 
on E^+ and Eg are different. As a result, by studying the energy spectrum of the final particles, it will be possible 

to extract + I-B2P and |Gi|^ + |G2|^ + 16(|Gip -I- IG2P). If in addition to the energy spectrum, the angular 

distribution of the electron relative to is measured, it will be possible to extract the parity violating combinations 
— |i?2p and I Gi p — I G2 1^-1-16(1 G2P — |Gip). Of course the larger the polarization of the initial muon, the 
higher the sensitivity of the angular distribution to these combinations. Thus, in principle by measuring the spectrum 
of the emitted particles and angular distribution of the final electron, it will be possible to derive the absolute 
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values of the couplings. The angular distribution of the final positrons also give information on |i3ip — and 
|Cip — IC2P + 16(16*2 P — 1^*1 P) (see the appendix). Thus, the following combinations can be measured from the 
study of the angular distribution plus energy spectrum measurements: 

\Bif, \B2f, \Cif + I6IG2I' and IC2I' + 16|Gi|2 . (7) 

However, without sensitivity to the spins of the final particles, it is not possible to derive information on the phases 
of the couplings. In sect. Ill A) we explore what can be learned from the polarization of the emitted electron. In 
sect. IIIBi we discuss the polarization of the positron. 



Polarization of the emitted electron 



Consider the decay fi'^ 



e^e e ' 



in the rest frame of the muon, where the final electron makes an angle of 9_ 



with the spin of the initial muon. Let us take the T„ direction parallel to the momentum of the emitted electron and 



(see Fig.[T]). The spin of the electron is determined by (ce- d^- ) with (|ce- P + |(ie- |^)^/^ = 1. 



That is the components of the spin of the electron are 

fg" • Se- = |ce- P - 14- • Sg- = 2Re[c*_4-] and ■ s^- = 2Im[c*_4-]. 

From the formulas in the appendix, we find that the differential decay rate in the rest frame of the muon is 



(8) 
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1287r4 
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J2 \M\^E^+Ee{mf, - - Ee)dE^^dE, 



+ (|Ci|V2 + 8|G2ni4-P + (|C2|V2 + 8|Gip)|c,-p] 
,-|V3+(|Gi|V2 + 8|G2p)|4-p-(|C2|V2 + 8|Gin|c, 



Re[G2G2*d:-Ce-]) 



(9) 



where is the polarization of the initial muon and dfl represents the differential solid angle of the momentum of 
the electron relative to the spin of the muon. To obtain this equation, we have summed over the spins of the emitted 
positrons and integrated over the azimuthal angle that the plane containing the momenta of these positrons makes 
with the plane of the spin of the anti-muon and the momentum of the electron. See the appendix for the details. 
Notice that there is no interference term between the chirality-flipping and chirality-conserving couplings. 
From Eq. in the appendix, we find that the longitudinal polarization of the electron is 
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(10) 



where 



p; 



Pi 
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P7 = 



2 

2 
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, \Ci? 



IG2I2 



8IG2P 



8|Gi|2- 
f 8|Gi|2 
■ 8|Gip -h a\i^2 

-8|Gip + 8IG2 



8IG2P 



|G2'^ 



(11) 



{sj,- ) is sensitive only to the absolute values of the couplings. Notice that by measuring {sj,- ) and its angular 
distribution, one can derive the same combinations that are listed in Eq. ([7]) i.e., the combinations that can be 
extracted from the angular distribution and energy spectrum (without measuring the spin of the final particles). 
Derivation of these combinations from the measurement of {sj.- ) can be considered as a cross-check for the derivation 

from the energy spectrum. Notice that even in the P^ = limit, (s^-) is still non- vanishing and provides us with 
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information on the parity violating combination Pi. That is while in this limit, the angular distribution of the electron 
is uniform and has no sensitivity to parity violation in the effective Lagrangian ([3]) (see Eq. (O). 
The transverse polarizations are 



and 



dT/dn\^^_^,,^, -dT/dn\ 

d_._=l/\/2 rf^_=-l/\/2 
ci^_=l/V2 !i^_=-l/\/2 

Re [(2.6)Bi*B2 - ^G\C^ - 4G2C2*] sin 6l_l 
P3- +P4-P^ COS0_ 



dT/dn\^^_^„^,+dT/dn\ 

Im [(2.6)Pi*P2 - 4GtCi - 4G2C2*]] sin^ 
+ P4"P^cos6'_ 



(12) 



(13) 



The transverse polarizations are maximal for electrons emitted in the direction perpendicular to the spin of the muon: 
i.e., 9- — 7r/2. Notice that (sj,-) and (sj,-) provide independent information on the real and imaginary parts of the 

same combination; i.e., (2.6P1P2 — 4GJCi — 4G2G2 ). Notice that both the moduli and the phase of this combination 
contain extra information beyond the combinations listed in Eq. ([7]). In the context of various models Gi and C'i 
can vanish. For example, as explained in the previous section, if the effective Lagrangian comes from integrating 
out neutral scalars, the Gi couplings will vanish. Within such models, (sj,-) oc RefPiPj] ''^^'^ {^t-) °^ Im[PiP2]- 
Considering that \Bi\ can be independently measured, the derivation of arg[PiP2] from both (sj,-) and (s^~) can be 
considered as a cross-check. 

In the limit of unpolarized muon (i.e., —>■ 0), (srp-) and (srp-) vanish. This is understandable because in the 

= limit, there is no preferred directions so we cannot define and r2^ (see Fig. [T]). Thus, to derive the 
CP-violating phase, arg[PiP2]i it is necessary to have a source of polarized muon which is quite feasible. For example 
if the muons are produced by the decay of pions at rest, they will be almost 100 % polarized. In fact, this is the case 
for the on-going MEG experiment. 

By measuring the polarization of the emitted electron, it is not possible to derive the values of all the CP-violating 
phases of the Lagrangian ([3]). Measurement of the angular distribution of the positrons does not provide any further 
information. As we shall see next, the transverse polarizations of the emitted positrons provide complementary 
information on the phases. 



B. Polarization of the positron 

Consider the decay /i^ — > e^^e^e^ in the rest frame of the muon where one of the final positrons makes an 
angle of 9^ with the spin of the initial muon. Let us take the direction parallel to the momentum of and 
f + EE f + X Sf,/\f+ X s^l (see Fig.[T]). The spin of e+ is determined by (c^+ d^+) with (|c^+p -t- M^+P)^/^ = 1. The 
components of the spin are 

■ ^et = Kt I' - Mef P. T+ ■ + = 2Re[c:+de+] and f+ ■ s^^ = 2Im[c:+d^+] . (14) 

111 1 "-1 1 1 '-'i 1 

From the formula in the appendix, we find that the differential decay rate in the rest frame of the muon is 

COS 6+ 



TTl 



3 

+ |PiP + \B2? + [|PiP(|c,+ P - M,+ |V3) - m\\d^^? - |c^+|V3)]P^cos0+ 
+ 16(|G2p|c^+p-H|Gind^+|2) 
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\Gi\'\d^ 



1 cos 6'+ 



+ 24(Rc[GiB*c*+d^+] -Rc[G2B2*c^+d*+])P^sin0+ 

where dfi represents the difTerential solid angle of the momentum of Ci relative to the spin of the muon. 
From the above equation, we find that the longitudinal polarization of is 
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where 
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(17) 



The transverse polarizations are 
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(18) 



and 
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; cos 6'-(_/3 



(19) 



Like the case of the electron discussed in sect. Ill Al the longitudinal polarization, {srp+), does not give information 
on the CP-violating phases and can be used only as a cross-check for the derivation of the combinations listed in 
Eq. d?]) by the methods discussed earher. The ratio of {Sj,+ ) and (sj.+ ) gives arg [B1C2 + B^C'i + 2AG\Bi - 24G2i?2]- 
Considering that the absolute value of this combination is also unknown, \{srp+)\'^ + \ {sj,+ )\'^ provides an independent 
piece of information. We have integrated over </>_ which means the measurement of the direction of the emitted 
electron is not necessary for this analysis. 



III. CONCLUSIONS AND PROSPECTS 



A large variety of the beyond standard models predict a sizeable rate for fX^ — > e"'"e~e+ exceeding the present 
experimental bound. In principle, by studying the energy spectrum of the final particles and their angular distribution, 
it is possible to derive the form of the terms in the effective Lagrangian leading to this process and extract information 
on the absolute values of the couplings (see Eq. ([7])). The effective Lagrangian responsible for /i+ e^e~e^ can 
include new CP-violating phases. In order to derive information on the CP-violating phases, we have suggested to 
measure the polarization of the emitted particles. In this letter, we have focused on the effective Lagrangian in Eq. ^ 
that can result from integrating out heavy scalar fields with LEV couplings at the tree level. The rest of the terms {i.e., 
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and Aji) are expected to be loop suppressed and are neglected in this analysis. We have shown that the transverse 
polarization of the emitted electron in /j,"*" — > e"'"e~e+ is sensitive to arg[2. 6-61-82 —iGlCi — 4G2C2] [see Eqs. (I12I13|) ]. 
That is while the transverse polarizations of the emitted positron is given by arg[i?iC2 + -BJC'i + 24C?Ji3i — 24G2-B2]- 
From Eqs. (|12ll3ll8ll9p . we observe that if the initial muon is unpolarized {i.e., — 0) the transverse polarizations 
of the emitted particles vanish. Thus, in order to derive the CP-violating phases, a source of polarized muons is 
required. 

In sum, the effective Lagrangian in Eq. ([3]) includes six new couplings and five physical phases. By measuring the 
energy spectrum of the final particles and the angular distributions relative to the spin of the initial muon, one can 
derive the CP-conserving combinations listed in Eq. ([7]): i.e., four out of the six CP-conserving quantities. Neglecting 
the loop suppressed and An couplings, the angular distribution cannot provide information on the phases. The 
transverse polarizations of the emitted particles provide four independent pieces of information on the phases and 
couplings. This information is not enough to reconstruct all the couplings but considerably reduces the degeneracy 
in the parameter space. 

We have also discussed the longitudinal polarization of the emitted particles. The longitudinal polarizations do 
not depend on the phases of the effective couplings. It is noteworthy that even in the = limit, the longitudinal 
polarizations of the positron, {srp+), and the electron, {sj,-), are nonzero and respectively yield information on the 

parity violating combinations jCip - jCsp + I6IG2P - 16|Gip and |Gi p - jGsp + IGlGsp - 16|Gi|2 + 2(|S2|2 - |Si p). 
Remember that in the P^ = limit, there is no preferred direction so the angular distribution of the final particles is 
uniform and does not yield information on the parity-violating combinations. 

There are running and/or under construction experiments that aim to probing signals for /i 67 [l^ and /i — e 
conversion on nuclei several orders of magnitudes below the present bounds on their rates. However, as shown 
in it is possible that while /i"*" e~^e~e'^ is round the corner, the rates of /i — > 67 and fi — e conversion are 
too low to be probed. In fact as shown in 4], the three experiments provide us with complementary information 
on the parameters of the effective LEV Lagrangian. If the muons are produced from the decay of pions at rest (like 
the case of the running MEG experiment ^12,]). the initial muons in — > eee will be polarized. On the other hand, 
there are well-established techniques to measure the polarization of the emitted particles in this energy range. In 
fact, such polarimetry has been used to measure the Michel parameters since 80's As a result, if the rate of 

/i+ — > e~^e~e'^ is close to the present bound and a hypothetical experiment finds statistically large number of such a 
process, performing the analysis proposed in this paper sounds possible. 

In this letter, we have focused on the LEV three-body decay of the anti-muon, /i+ e+e^e+. The same discussion 
applies to the decay of the muon, /i^ —>■ e^e+e^. In this mode, the transverse polarizations of the electrons would 
give arg[-BiG| -I- -BJCi -I- 2AG\Bi — 24G2-B2] the transverse polarizations of the emitted positrons would give 
arg[2.6i?i-B2 — 4GiGi — AG2C2]. The method of measuring the polarization described in 14] is based on studying 
the distribution of the photon pair from the annihilation of the emitted positron on an electron in a target. If this 
method is to be employed, only the polarization of the positron can be measured. Thus, to derive both combinations, 
the experiment has to run in both muon and anti-muon modes. 

The three-body LEV decay modes of the tau lepton such as t ^ eee or t — > can also shed light on the 
underlying theory. Recently it has been shown that by studying the angular distribution of the final particles in 
T — > MMM 3't the LHC, one can discriminate between various models [l5| . Discussions in the present letter also apply 
to the decay modes t — > eee and r — > /i/l/i. 



Appendix 

In this appendix, we derive the decay rate of /i^ — > e^e~ej. We first derive the decay rate into an electron of 
definite spin, summing over the spins of e^ and e^. We then concentrate on the spin of ef and derive the decay rate 
into a positron of definite spin, summing over the spins of the electron and the other positron. 

With effective Lagrangian ([3]), we find 

+ + - u -1±75 - 1 - 75 p -1+75 - 1-75 

M{^r ^ e{e ej) ^ Bi^l — ^ — eie — ^ — e2 - Si/i — ^ — e2e — ^ — ei 

, p-l-75 -1 + 75 j,-l-75 -1 + 75 
+ S2/i — ^ — eie — ^ — e2 - -82/1 — ^ — e2e — ^ — ei 

^-1 + 75 -1 +75 ^-1 + 75 -1 + 75 

+ '-'iM — 2 — — 2 — — ^ — e2e — ^ — ei 
, ^-1-75 -I-75 ^-1-75 -I-75 
+ G2/i — 2 — ^1*^ — 2 — ~ — 2 — — 2 — 

.^-1^75o*t1 + 75 .^-1+75o*t1-75 

- 4Gi/ic — - — 7''e e^ c — - — e2 - 4G2//C — - — 7"e e^ c — - — e2. (20) 
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Decay rate into e with a given spin 



Consider the decay /i+ e^e^e^ in the rest frame of the muon. Since we are interested in the spin of the 
electron, it is convenient to use the coordinate system defined as = P^- /\P^-\, T2 = (2^3^ x Sfj_)/\Tr^ x and 
= X Tr^ . In this coordinate system, 

P,,+ = (m^, 0,0,0) P,+ = (1, sin a cos sin a sin(/)_|-, cos a) 

Pe- = (Pe,0,0,Pe) P(.+ = {'>n^ — E^+ — Pg , — sin a COS (/)+ , — £'g+ sin Q! sin 0+ , — i?e — COS a) (21) 

where the electron mass is neglected (see Fig.[T]-a). Writing the kinematics and neglecting effects of 0(mg/m^) <^ 1, 
we find 



cos a 



2E + Ef. 

Summing over the spins of the emitted positrons, we find that 

r.27r 



(22) 



(27r)4 / (2E^^){2EA ^ |Afpd^+ 



1 2 



= (iPinc,-p + iP2pK-ni?e 



(iPinc,-|2-iP2ri4-i')i?, 



(m^ - E^+{\ - cosa)) + (1 - cosQ;)(mp - Pe - E^-v) 
cosq![to^ — £'g+(l — cos a)] — (1 — cosQ;)(i?g+ cos a + E^) cos0_ 
+ [(|Ci|2 + 16|G2|")|de^|'(l + P^cos(?_) + (IC2I' + 16|Gi|2)|ce-|'(l -Ppcos0_)] P,+ K - Se(l -cosa)] 
- 2P^Re[BiP^4_c*-]£:^+(l-cosa)(P^+(l-cosa)-m^)sin6'_ 

+ 8Pp(Re([GiCi*4-c:-]+Re[G2C2*d:-cl-])P,+ (Pe(l-cosa)-m^)sin0_ , (23) 

where (cg- d^-) determines the spin of the emitted electron [see Eq. ([8|)]. 

The differential rate of the decay into an electron in a direction that makes an angle of Q- with the spin of the 
initial muon is 



E 

s + ,s + 

=1 "1 



dT{iJ,+ ej ^e-e+^ 



Jo JraJ1-E,^^^^^ 



|M|2p,S^+(m^ - Pe - E^+)dE^+dE,dq^+, 



where dVl is the differential solid angle determining the orientation of the emitted electron. The factor Ei,E^+{inf^ 
Ee — E^+) comes from the momentum-space volume for a three body decay [i.e., from i5*(P/i — Pe — P^+ 
P^+)d^Ped^P^+d^P^+]. Inserting \M\'^ from Eq. we obtain 



E 



(ir(/i+ e{e 63 ) 

dn 



^1 ^1 



8(27r)4 

p^cos0_(iBinc,-|2-|P2n4-p; 



1/2 /.1/2 



Jl/2-y 
1/2 <.l/2 



-1 + + 4a; - 2y^ - 4a;^ - 4:Xy)dx dy 
2x + 2y-l 



[l-2x-2y + 2xy 



Jl/2-y 



y 



+ ((|Gi|2 + 16|G2|')|4-|'(l+PpCos0_) + (|G2|' + 16|Gi|2)|ce-|'(l-P^cos(?_)) / / y{l-2y)dxdy 

Jo J1/2-V 

/•1/2 rl/2 

- 8P,,(Re[GiGi4-c*-] +Re[G2G2*d*-Ce-])sin6l_ / / y(l - 2y)d 

Jq J1/2-V 



{l-2x-2y + 2xy + y^)]dx dy 
1/2 .1/2 

1/2-y 



ix dy 



f'^ f-'^ 2x-l 
- 2Pf,Rc[BiB;de-cl-]sm9- / (2a; + 2y - 1) dx dy 

Jo Jl/2~y y 

(24) 

The integrals are all finite and lead to the numbers in Eq. ^ 
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Decay rate into ef with a given spin 

Let us now concentrate on one of the positrons, e^, whose momentum makes an angle of 6^ with (see Fig. 
[T]-b). To perform this analysis, it is convenient to work in the following coordinate system: = /\P^+\, 
T2 = X s^/\T^ X Sfj\ and = T2 x . In the rest frame of the muon and in this coordinate system, 

=(to^, 0,0,0) Pp- = i?e(l,sinQ!Cos0_,sinQ!sin0_,cosa) 
P.+ = 0, E^+) Pg+ = (to^ - ~ Ee, ~Ee smacos(p^,~Ee smasiiKj)^, -E^+ - E^ cos a) (25) 



where a is given by Eq. (p2|) . 

Summing over the spins of and e~, we find 

{2nf{2E,){2E^^) J ^ |Mpd0_ 



S _ ,S -I- 

|2| J |2 



= i\Bi\'K+\' + \B2\'\d^^\')E4m^ - E^+il-cosa)] + m\'\d^^\' + 
+ 16(|G2p|c^+P + \Gi\^\d^^\^)Ee{mf, - E,{\ ~ cosa)) 
+ {\Ci?K+? + \C2?\d^^\hEe{m,, - E,{1 - cos a)) 

(|Bl|2|c^+|2 _ |B2pMe+nK - ^e+(l - COS a)] COS 0+ 

+ P^dPaPlc^+P - |Bind^+P)(l -cosa)(£;eCosa + i;,+ )l EeCose+ 



|S2p|c^+P)(l -COSa)Se(m^ -Pe 



1^) COS sina(l - cos a) sin^ 



+ P^dCiHc^+p - IC2 nd^+n£;e cos a(m^ -i;e(l- cos a)) cos 0+ 
+ leP^dGaPlc^+P - \Gi\^\d^+\^)Ee{m^ - i;e(l - cosa)) cosacos6l+ 
+ PpRe[-BiC2Cg+d*+]m^-Be sin6'+(l + cos a) 
+ P^,Re[B2C*c*+d +]mfj,Ee sine'+(l + cos a) 

Re[GiBlc\d^+] -Re[G2B;d*+c^+] \ Ee{m^, - Ee{l - cos a)){l - cos a) sin 9+ 



(26) 



where {c^+ d^+) determines the spin of the emitted positron [sec Eq. p^ ]. The differential rate of the decay into a 
positron in a direction that makes an angle of 0+ with the spin of the initial muon is 



E 



d^(^^ 



dn 



2tv /•ni^/2 
^0 



J2 \M\^E,E^+ (m^ -E,- E^+)dEe dE^+ 



where dfJ is the differential solid angle determining the orientation of . Inserting \M\^ from Eq. (j23p . we obtain 
(ir(/x+ — > e~ej) 



E 

8(27r)4 



dr2 

(Me+n^ll' + IVl'l^^P) 

1/2 <.l/2 



1/2 ^.1/2 



1/2-y 



(1 - X - ?/)(2x + 2y - l)da; 



+ c 



,,\W + K^\W) / 2;(l-2y)da;dy 



+ i\c\f\c^ff + \c2fKfn 



i/2-y 

1/2 (.1/2 



'1/2-y 
"1/2 fl/2 



x{l ~ 2x)dx dy 



+ 16(|G2nc^+P + |Gind^+|2) / / x{l-2x)dxdy 

^ ^ Jo Jl/2-v 

+ p^cos0+(isind^+|2-|i32nc,+|2) 



l/2-iy 

1/2 ^.1/2 



1/2-y 

1/2 (.1/2 



1 — 2a; — 2-!/ + 2a;?; , , , 
(2a; + 2y- \){-y ^ -)]dx dy 



.cose+(|Cinc,+ p-|G2nd,+ n / / (l-2a;) 

^ ^ Jo Jl/2-v 



1 - 2a; - 2y + 2xy 
2^ 



2?/ 
da; 
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\cose+{\Bi\'\c^^\' -\B2\^\d^^\') [ ^ [' y{l-2y)dxdy 

^ ^ Jo Jl/2-v 



ll/2-y 



2P^(Rc[SiC*c +Rc[B2C*c%d +])sin6l+ / / 

^ Jo Ji 



Jl/2-y 2j/ 
1/2 1-1/2 

(2xy + 0.5 — X — y)dx dy 

/2-y 



4P^(Rc[GiS*c*+d^+] -Rc[G2B2*c^+d*+])sin0+ I ' I ' 2x){2x + 2y I) 

^ ^ ^1 Jo •ll/2-y 4X 



(27) 
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